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ABSTRACT 

^ ! The next generation of neutrino and 7-ray detectors should provide new in- 

O . sights into the creation and propagation of high-energy protons within galaxy 

^ ' clusters, probing both the particle physics of cosmic rays interacting with the 

background medium and the mechanisms for high-energy particle production 
within the cluster. In this paper we examine the possible detection of 7-rays (via 
the GLAST satellite) and neutrinos (via the ICECUBE and Auger experiments) 
^ ! from the Coma cluster of galaxies, as well as for the 7-ray bright clusters Abell 

85, 1758, and 1914. These three were selected from their possible association with 
unidentified EGRET sources, so it is not yet entirely certain that their 7-rays are 



O 
H 



c/2 ! indeed produced diffusively within the intracluster medium, as opposed to AGNs. 

It is not obvious why these inconspicuous Abell-clusters should be the first to be 
seen in 7-rays, but a possible reason is that all of them show direct evidence of 
^ [ recent or ongoing mergers. Their identification with the EGRET 7-ray sources 

Q\ [ is also supported by the close correlation between their radio and (purported) 

^ i 7-i'ay fluxes. Under favorable conditions (including a proton spectral index of 

I 2.5 in the case of Abell 85, and ~ 2.3 for Coma, and Abell 1758 and 1914), we 

O I expect ICECUBE to make as many as 0.3 neutrino detections per year from the 

Q I Coma cluster of galaxies, and as many as a few per year from the Abell clusters 

85, 1758, and 1914. Also, Auger may detect as many as 2 events per decade at 
I ~ 10^^ eV from these gamma-ray bright clusters. 

Subject headings: acceleration of particles — galaxies: clusters: individual (Coma) 
— neutrinos — radiation mechanisms: non-thermal — relativity — X-rays: 
galaxies 
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1. Introduction 

A population of high-energy cosmic-ray (CR) protons, pervading galaxy clusters, was 
first invoked (Dennisson f 980) to explain highly polarized radio luminosities {Lr ~ f 0*^*^ — f 0^^ 
erg s~^) now observed in over 30 clusters. Evidently the radio light is synchrotron emission 
from highly relativistic electrons (7e ~ 10^) gyrating in large magnetic fields (S ~ 1 /xG), 
despite the fact that such electrons would lose their energy to radiation well before they 
could traverse the ~Mpc extent of observed emission. However, protons injected into the 
cluster diffuse throughout the intracluster medium, essentially unchanged and confined for 
cosmological times. In this environment they may interact with the background gas and 
create a decay cascade. Among the products of this decay are the relativistic electrons 
responsible for radio emission. 

Since Dennisson's (1980) pioneering work, conclusive observations of this 'secondary' 
process for electron production within clusters have been few. In Hercules A (Nulsen et al. 
2005b), MS 0735 (McNamara et al. 2005), and Hydra A (Nulsen et al. 2005b), vacated areas 
of X-ray bremsstrahlung indicate the background of thermal electrons is being pushed aside 
by relativistic electrons injected from the axes of a central source, implying a phenomenally 
large power, Lp ~ 10^^ erg s~^. These 'bubbles' coincide with sources of radio emission 
as the protons produce a cascade of synchrotron-emitting electrons, confirming that the 
secondary decay process does occur in clusters. The protons released during such an episode 
would be confined and built up within the cluster over cosmic times, totaling some 10^^ 
erg of nonthermal hadrons — very near energy equipartition with the thermal background 
(Berezinsky, Blasi, and Ptuskin 1997; Hinton & Domainko 2006). This is, we find, the order of 
nonthermal energy required to be deposited in the nonthermal hadron population for proton- 
proton scatterings to produce 7r° decay 7-rays within the EGRET detection sensitivity. 

The energy liberated by the infall of the galaxy group NGC 4839 toward the center 
of the Coma cluster, observed by XMM-Newton (Neumann et al. 2001), is an example of 
another potential source of high-energy protons. NGC 4839 achieves a velocity of ~ 1,400 
km s^\ and since the sound speed corresponding to Coma's gas temperature of ~ 8 keV is 
~ 1,000 km s~^, the subcluster's supersonic motion is expected to produce shocks (which 
Neumann et al. claim to observe directly in the imaging of NGC 4839). It is hkely that 
the shock compression in these regions energizes a fraction of the thermal particles within 
the intracluster medium (ICM) by first-order Fermi acceleration. However, how the protons 
are actually energized is not directly relevant to our analysis. For example, since the cosmic 
rays remain trapped within the cluster for over a Hubble time, it is likely that second-order 
acceleration by a turbulent distribution of Alfven waves may be even more effective at helping 
the protons attain their highest energies (see, e.g., Liu, Petrosian, and Melia 2004; Liu, Melia, 
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and Petrosian 2006). Our principal goal is to determine the observational signatures such 
a population would produce, should it be present within the cluster. But cluster mergers 
do provide us with a viable mechanism for energizing the cosmic rays, regardless of how the 
individual particles are ultimately accelerated. In this regard, we note that the gravitational 
energy ~ 10*''' erg released in a large relic merger could reasonably supply ~ 10^^ erg in a 
nonthermal hadron population, again sufficient to supply detectable 7-rays. In this case, the 
radio halo may be interpreted as a signature of the cosmological development of a cluster 
(the oldest known cluster being some 9 billion years old) . 

Aside from a radio halo produced by secondary electrons, conclusive evidence for the 
properties of CR protons within galaxy clusters rests on the detection of other products of 
the protons' decay cascade: 7-rays and neutrinos. The closest active galaxy cluster, the 
Coma cluster, does not emit copiously enough to have been detected by EGRET given 
its ~ 10"^^ cm~^ s"^ flux sensitivity, although we will show that this does not preclude a 
possible detection by the GLAST satellite or a detectable neutrino flux. Our preliminary 
identification of clusters which are 7-ray bright, and, therefore, plausible extrasolar neutrino 
candidates, has come from a systematic search for likely clusters among unidentified EGRET 
sources. 

In a search for new radio galaxy clusters, the NRAO VLA (Giovannini et al. 1999) 
and Wcstcrbork Northern (Kcmpner & Sarazin 2001) sky surveys were systematically cross 
referenced with X-ray bright objects from the ROSAT survey (Ebcling et al. 1998), providing 
a list of more than 30 known cluster halos or relics, presented in Table 1. While none of 
the clusters EGRET was hoped to detect produce a measurable 7-ray fiux (Reimer et al. 
2003), three radio sources coincide with unidentified 7-ray bright objects detected by EGRET 
(Colafrancesco 2002)— the cluster Abell 85 and the EGRET source 3EG J0038-0949, Abell 
1914 and EGRET source J1424+3734, and Abell 1758 and EGRET source 3EG J1337 +502a. 
These three EGRET sources are not unequivocally the quoted clusters, but here we shall 
assume they are in fact the same objects. 

The assumption is not certain. In a study by Reimer et al. (2003) it was shown that 
the probabihty that one of the 170 unidentified EGRET sources and a member of their list 

of 58 considered galaxy clusters coincided was 48.1%. Of their list of 58 clusters worthy of 
consideration, only Abell 85 and the unidentified source J0038-0949 showed a considerable 
overlap, which could therefore be explained as an entirely chance occurrence. However only 
six of their 58 considered clusters — included for their proximity {z < 0.14) and X-ray 
brightness — were shown by the NRAO VLA or Westerbork Northern radio sky surveys to 
be radio halo/relic clusters, so that in only these six cases is there any reason to believe a 
significant nonthermal proton population may be present. Also, Abell 1758 and 1914 were 
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not considered by Reimer et al. (2003), as they are above z > 0.14. 

The three candidates are not the most X-ray bright clusters, most of which were assigned 
only upper limits in the nine years of EGRET observations. It is not obvious why three 
inconspicuous Abell-clusters are the first to be seen in gamma rays when the closest and 
most well-studied clusters — known to have evidence for non-thermal particle populations 
from radio, EUV, or hard X-ray detections — are not. Perhaps a key is that each of the 
7-ray bright candidates shows direct evidence of recent or ongoing mergers, indicating (if 
confirmed) that this may outcompete the AGN outburst as a source of cluster cosmic rays. 
In Abell 85, a 4 Mpc filament (Durret et al. 2005) forms a chain of several groups of galaxies. 
In particular, an impact region (the 'south blob') displays a temperature (8.8 keV) far in 
excess of the rest of the cluster (2.7 keV). In Abell 1914 (Govoni et al. 2004) a NE-SW 
arch-like hot region extends across the cluster center and may represent an ongoing shock. 
Abell 1758 is a double cluster (David & Kempner 2004): A1758N is in the late stages of a 
large impact parameter merger between two 7 keV clusters, while A1758S is undergoing a 
smaller impact parameter merger between two 5 keV cores. In each case, the fiux of 7-rays 
correlates with the radio fiux, as it would if both were linked to secondary decay products. 
Abell 85 also shows evidence of a hard X-ray excess (Lima-Neto et al. 2001), indicative of a 
nonthermal population of electrons. 

In the end, however, we caution that while the identification of Abell clusters 85, 1758, 
and 1914 with three unidentified EGRET sources may be motivated by the reasons we have 
outlined here, it is not entirely clear that their 7-rays are indeed produced diffusively within 
the IGM. Thus, one must accept the results presented in this paper, particularly the predicted 
neutrino fluxes, with this important caveat in mind. 

As clusters confine CR protons over cosmological times, they provide a promising source 
for a possible first detection of high-energy (> TeV) extra-solar neutrinos. Aside from 
probing high or ultra-high energy particle physics processes occurring within astrophysical 
sources, any such neutrino (or 7-ray) detection would also establish the properties of the 
proton population within the cluster and settle the possible role of such populations in various 
cluster phenomena — e.g., generating the radio halo (Dennisson 1980), their contribution to 
the 'cooling fiow' problem (Silk 1995), or their acceleration within cluster merger shocks 
(Sarazin 2004). 

In this paper we take the optimistic view that protons within the Goma cluster create a 
flux of 7-rays which is within the detection hmits of the GLAST satellite, without violating 
the EGRET null detection. In the same spirit of optimism, we identify the Abell clusters 
85, 1758, and 1914 with their gamma-ray bright counterparts. The resulting predictions 
for neutrino and gamma-ray detections for these four clusters will determine whether this 
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optimism was well-founded. 

Work on the hadronic origin of one or more spectral components in the overall clus- 
ter emissivity has been quite extensive, and our investigation overlaps with several earlier 
treatments, though our model is unique in the treatment of the particle kinetics and proton- 
electron scattering events. Dolag and Ensslin (2000) considered the possible origin of the 
radio halos from a hadronic secondary electron injection, based on a detailed modeling of 
the cluster magnetic field (see also Blasi & Colafrancesco 2000, Jones 2004, Brunctti ct al. 
2007, Marchcgiani et al. 2007). Subsequently, Pfrommer and Ensslin (2004) constrained the 
cosmic-ray population in several cooling fiow clusters from their radio and 7-ray cmissivities 
(see also Reimer et al. 2004, and Sarazin 2007). Not surprisingly, the 7-ray emission from 
proton-proton collisions in clusters may contribute a non-negligible fraction of the extra- 
galactic background radiation field (Kuo et al. 2005). And though we will not be addressing 
7-ray line emission in this paper, cosmic rays may also collide with heavier nuclei in the ICM 
to produce measurable 7-ray line intensities (lyudin et al. 2004). An extensive review of the 
7-ray emissivity in clusters may be found in Blasi et al. (2007) . 

The issue of hadronic acceleration during large-scale structure formation has been con- 
sidered by several authors, including Miniati et al. (2001), Berrington & Dermer (2003), 
Blasi (2004), Ensslin et al. (2007), Jubelgas et al. (2008), Ando and Nagai (2008), Nakar ct 
al. (2008) and, in the context of secondary acceleration due to stochastic processes, Brunctti 
ct al. (2004), Brunctti and Blasi (2005), and Pctrosian and Bykov (2008). But cosmic 
rays may also play an important dynamical role throughout the cluster's lifetime. Collisions 
between energetic hadrons and the ICM have been invoked as a means of heating the gas 
condensing toward the middle of the cluster, thereby preventing an overly rapid cooling of 
the fiow (Colafrancesco, Dar, and De Rujula 2004, Guo & Oh 2008). 

Our goal in this paper is specifically to calculate the neutrino fiux expected from 7-ray 
bright clusters, but our approach differs from earlier attempts in several important ways, 
which we will describe over the next several sections. Our model requires that the number of 
protons be consistent with a secondary explanation of the Coma cluster's radio halo, using 
magnetic fields (~ l/^G) consistent with observations of Faraday rotation, rather than those 
(~ 0.1/uG) consistent with an inverse- Compton description of the hard X-ray excess detected 
by Beppo5'^X and RXTE (see also Wolfe and Meha 2006). The only model-dependent 
characteristics of this work are then the distribution of gas and that of protons diffusing 
throughout the cluster, which — together with the details of calculating the various links 
in the decay cascade — we discuss in Section 2. Our technique for calculating the detected 
number of neutrinos, given the predicted fiux of Coma, is discussed in Section 3. Finally in 
Section 4, we review the forecast for those clusters for which a positive 7-ray detection has 
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been proposed. 

2. Model Chciracteristics and Flux Calculation 

2.1. Emission spectra of nonthermal plasmas 

Our model first requires the distribution of electrons, neutrinos, and 7-rays produced 
by proton-proton collisions and subsequent decay of secondaries. Gamma rays are produced 
from the decay of neutral pions pp ^ ti^ ^ 27 at a characteristic energy of ~ 70 MeV in the 
pion rest frame; a similar number of neutrinos arc also produced via — tt^ — > /i^ + z/^(z/^), 
and in the subsequent decay of muons into electrons and positrons, /x^ — > e"^ + +i'e{ue) 
(given typical power-law parent proton distributions, neutrinos from neutron /3-decay may 
be neglected in these calculations). Both neutrinos and 7-rays may also be produced via 
proton interactions with background light fields. For the cosmic microwave background, pion 
production is characterized by a peak at 6 x 10^^ eV — approximately the energy of a 140 
km/h tennis ball — while for infrared light emitted by galaxies and proposed (Blasi 2005) as 
a dominant source, this peak shifts to 2.5 x 10^*^ eV. All three decay sources probe different 
energy regimes: electrons, via synchrotron radio (10^^ eV) or HXR inverse- Compton (10 eV); 
7-rays, at 70 MeV and notably detectable by GLAST and EGRET; and TeV+ neutrinos, 
which directly probe the tail of the proton distribution (10^^-10^^ eV) via pp scatterings as 
well as its extreme limit (> 10-*^^ eV) via interactions with the CMB. 

Our recipe is to assume or imply (via radio observations, see Sec 3.2) a proton distri- 
bution, which in a practical sense means setting a spectral index and normalization to a 
power-law np{Ep) = UpoE'p' . From this alone follow two injections of pions — one for proton 
interactions (pp) with hydrogen in the ICM, and one for interactions with background light 
(p7). The observationally relevant decay products of these pions are electrons, 7-rays, and 
neutrinos (see also Blasi and Colafrancesco 1998), which for the particles of interest may be 
found by summing over the kinctically allowed energies in each link of a decay chain; this 
section is devoted to a description of this process. 

For proton-proton interactions, the neutral pion emissivity is 

ql',{E,o) = CUH [ dEpUpiEp)^-^^^^^ , (1) 

and similarly for charged pions, via an altered cross section. The threshold proton kinetic 
energy for production of a pion in a proton- proton collision is Ethres = 2m.,^{l + rriT^/Amp). 

Here we take the cross-section for pion production for energies below the isobar reso- 
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nance {Ep — 7 GeV) from Dermer (1986) — where also each of the various decay channels 
is treated — while above the resonance we take the cross-section to be a simple scaling form 
used by Blasi & Colafrancesco (1999), and Fatuzzo & Melia (2003). We note in this context 
that while collider data at these extreme energies do not exist, the proton-proton collision 
cross-section has been examined by the Fly's Eye detector, which implies from atmospheric 
fluorescence that the collision cross-section violates the scaling assumption when ^/s is tens 
of TeV. Unfortunately the pion production cross-section is determined both by the rate of 
proton collisions, Rpp — uncdapp/dEp, and by the varying multiplicity, M{Et^, Ep), of pion 
production with proton energy (Moffeit et al. 1972, Markoff and Meha 1997, 1999). For 
this we again depend on accelerator data, which exist only up to 1,800 GeV — so that the 
limitation on our knowledge of the cross-section is compounded by a lack of certainty regard- 
ing pion production. Simply, predictions and observations of neutrinos or 7-rays in clusters 
probe energy regimes not previously accessible. 

Once the pion source function is found, the neutrino, 7-ray, and electron fluxes follow 
from kinematic concerns discussed fully in Stecker (1979) for 7-ray decays and Marscher et 
al. (1980), and Zatsepin & Kuz'min (1962) for neutrinos. For decay photons we have 

r (i''''(E ') 

<^^-^-'l/'" iEl!-^^^n ' (2) 

where the minimum kinetically allowed pion energy is (p = E^ + m^oC^ / {AE^). Well above 
the 70 MeV region, the pion spectrum shares the proton parents' spectral index, s^r = Sp. 
For electrons, the source function is 



'^l - Je, dEe Je^^u /3„E^ JEthr.,{E^) dE^ 

(3) 

where dP/dE^. is the three-body decay probability (Fatuzzo & Meha 2003). Here, the hmits 
of kinematically allowed pion energies are 

IE 

Tpmin,max M f a\ 

- {l±(3,)-rml/ml{lTP,) 

(with minimum corresponding to the upper sign, maximum to the lower). The threshold 
Ethresh{.ET,)ioi produciug a pion with energy Et^ is similar to that in Eq. (1). 

For neutrinos, we must distinguish between particles produced by pion decay, and those 
produced by muon decay, since these are two- and three-body decays, respectively. For pions 
the source function is 

/■oo 

C(^-J = / ni.)QMdl. . (5) 

' m.i.n. 
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In the case of electron injection resulting from proton collisions we integrated over the three- 
body decay probability dP/dEe- Here /(77r) plays that role, as the spectrum of two-body 
decay 

Here, E^ ^ 29.8 MeV is the CMS energy of the emitted neutrino. In this instance the 
minimum kinetically allowed energy is 

7L„ = EI/2E, + E^/2El . (7) 
In the decay of charged muons, the muon neutrino flux is 

q>:;{E,) r f{E,) f g.(7.)|^ dx dj^ , (8) 

where electron neutrinos are produced at approximately half the rate of these muon neutri- 
nos. The interaction angle cosine, %, is kinematically limited to 

Xmm = -1 , for 3.68 < 3.557^ 

= ^13.54 - 12.672 , for 3.68 > 3.557^ , (9) 

while the lower limit for 7^ is 

1„ ^ AE, ■ ^ ^ 

In this case we require f{E^), the spectrum of three-body decay, which is 

f{E^) = 167^(1- 1(3 + /3J)C)C'— , forO<C<[l-W2 



5 4 /8C 



- 3(1 + (3^)) CM,m, , for [1 - /3J/2 < C < [1 + /3;.]/2 , (11) 



where C, — {t^^) and (3^ — f^/c. Since detection lies far above 70 MeV, it suffices to 
approximate neutrinos from pion decay as the same in number as those from subsequent 
muon decay. 

In the case of proton-photon interactions, we may use the expressions (2, 5, &8) unal- 
tered, changing only the source function of pion production. This is now given by (Bottcher 
& Dermer 1999) 

c [°° d'^a 

(iA17t) = ^J^ d-fpTipi-fp) J ^dx den^{e){l- (3px)j^-^, (12) 
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where Up and rij are, respectively, the proton and photon distributions. The result is a flat 
neutrino emissivity dN^/dE dV dt until threshold, with a spectral index one higher than 
that of the proton distribution afterwards, s^, — Sp + 1. Again, x is the interaction angle 
cosine. 

The cross-section for charged photopion production (Mucke et al. 2000, and references 
cited therein) features two prominent decay resonances (Fig. 1): the first, at — 330 MeV, 
has a height of 2.45 x 10~^^ cm^ and a decay width of 0.15 GeV. The second resonance, at 710 
MeV, has a decay width of 0.18 GeV — but, due to the steepness of the proton distribution, 
it docs not play a major role in the neutrino flux. Note that this height is some factor of 
1.5 greater than that assumed by Dermer and Bottcher (2003). The cross-section for neutral 
photopion production shows only a single resonance of height 3 x 10~^^ cm^ and a decay 
width of 0.15 GeV. Our fltting function for a resonance at mass M, of height (To and of width 
r is of Breit-Wigner form, 

^(^') = (e'-M2)2 + rV ' ^^^^ 

where e' = 7pe(l — /i) is the photon energy in the proton rest frame. We have tested the 
calculation with a simple 5-function resonance, however, and the results are nearly identical. 

Finally, we require the synchrotron emission from injected electrons to set the level of 
ambient relativistic protons. In practice this is evaluated numerically, but a sufficient guide 
may be given by the approximation for the electron injection (Mannheim & Schlickeiser 
1994), 

rrin f E^ \ 13 ( nip yeo-i _3 _^ 

= 707;^^-n70I^j^l2^-^^^"-^^^H24^^ ' 

= /r,„,-7e"'^° ■ (14) 

Assuming the dominant loss mechanism is either inverse Compton scattering with a black- 
body photon background {Tqmb — 2.73 K), or radio synchrotron, radiative losses are given 
by —dEe/dt = asE^, with the constant = {4/3)aT cne{ecMB + ^B)/fT^e(?- The energy 
density {ecMB) in the CMB dominates over that (e^) in the magnetic field. The equilibrium 
distribution of electrons due to injection against these losses is (Blumenthal and Gould 1970) 

n{le) = ^^ ■ (15) 

meC^as{Seo - 1) 

The radio synchrotron emissivity (in units of energy per unit volume, per unit time, per unit 
frequency) associated with this distribution is then 

^ l-157r^ _ o„ . TT — > (16) 



dV dvdt nieC^ as{seo — l)\uJ 
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where ub is the gyrofrequency. For a final handle on the electron population in clusters based 
on the presence or absence of a hard X-ray excess, the corresponding Compton scattering 
emissivity off the CMB (in units of photon number per unit volume, per unit time, per unit 
energy) is 

"^^^ = 1.8^ ,, (fercMB)^^-°+^^/^6-(--°+^)/^ , (17) 



dV de dt h^c^ mec'^as{seo — 1) 
where tq is the classical electron radius, and we have used the fact that the hard X-radiation 
is produced below the Klein-Nishina region to simplify the cross section. 



2.2. Cluster modeling 

We now have spectra or distributions for each of the links in the decay cascade — 
electrons, 7-rays, and neutrinos — and have derived expressions relating electron number 
to observed radio or nonthermal X-ray emission. It remains only to assume a radial distribu- 
tion for the background gases and for energetic protons diffusing throughout the cluster. For 
simplicity we will adopt a single point source of protons at the center of the cluster: realistic 
morphology, time dependence, or non-Kolmogorov proton diffusion through the cluster may 
seem to increase the accuracy of the model, but without a motivating observation only serve 
to make it less reproducible. Below we describe in detail how we calculate the 7-ray and 
neutrino fluxes from Coma, the results of which are presented graphically in Fig. 2. 

We take the background gas within Coma to be described by a /3 model (CoUess & 
Dunn 1996) with ricore = 3 x 10"^ cm~^, (5 ^ 0.7, and Vcore = 0.25 Mpc, for which 

r / r \ 2l -3/3/2 

n(r) = nc 1 + f — J . (18) 

While the thermal and nonthermal electron populations are well determined by X-ray and 
radio observations (Fig. 3), respectively, no similar tool directly determines the nonthermal 
proton population. Here, we assume high-energy proton injection takes place at a central 
source and these then diffuse outward according to (Blasi & Colafrancesco 1999) 

^"""^g^'^'^^ = D(E,) n,(E„ r, t) + Q(E,)5(r) + ^KE,)n,(E„ r, t) . (19) 

The highest possible energy attainable by protons is actually given (Dermer & Berrington 
2002) by the duration of acceleration, rather than any loss term b{Ep). Loss, and an addi- 
tional term representing proton escape from the cluster, we neglect. Assuming the cluster is 
relaxed, we then may integrate to give 
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where Qp{Ep) — QoEp and Qq is normalized to synchrotron observations. 

The resonant diffusion of protons in magnetic field fiuctuations with a power spectrum 
P{k) in wave number k is given by (Blasi & Colafrancesco 1999) 

where tlIp) — pc/eB is the Larmor radius. In this paper we take fiuctuations of the magnetic 
field to have a Kolmogorov spectrum, 

/ k \ -5/3 

m = ^(^o)(^) , (22) 
and assume the spectrum is normalized to 

P{k)dk ^ B^. (23) 



/ 



' ko 

Thus 

P(to) = -^^B' , (24) 

and 

D{p) = l/3cA;o"'^'(eS)-V3£;V3 _ (25) 

A typical value for the smallest scale on which the magnetic field is homogeneous, do — 1/ko, 
is > 20 kpc, since the magnetic field is being stirred by the virial motion of galaxies within 
the cluster, for which do = {^n/{3Ngai)Y/^Rciuster- 



3. Detection Counts 

In our model, the number of photons above detected with GLAST, writing L43 as 
the luminosity in units of 10^^ ergs s~^, is (Berrington & Dermer 2005) 

A^(> E^) ~ 35 V(^7/GeV)-^-°^ ^ , (26) 

where we use a cluster evolution time tyr = 0.7 Gyr, and L43 = 0.1 through 0.7 coinciding 
with Sp = 2.1 through 2.5. The detection significance rja for Coma, assuming it is a point- like 
source (i.e., the point spread function for a given photon energy is larger than the 18' angular 
extent of Coma) is (Berrington & Dermer 2005) 

77.~5.4Vv(^)(VGeV)°^ (27) 
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This yields a rja- = 5(7 detection significance. 

We may also easily construct an argument for why neutrinos may be detectable. Suppose 
the 7-ray flux detected from the 7-ray bright Abell clusters is around 

$^(1 TeV) ^ 10"" photons cm^^ (28) 

at a few TeV (see Fig 4). This implies a similar number of muon neutrinos. Following 
Gaisser & Stanev (1984), the flux of muons produced in deep-inelastic interactions between 
neutrinos and ice is nearly 

d^, ~ ^ 10-2° GeV-i cm-2 s'^ , (29) 



with the deep-inelastic cross section a, muon energy loss rate a, and Avogadro number Nav 
The muons leave the interaction with about l/3rd the incoming neutrino energy, so that 

^,^E,^^3xlO-''cm-'s-\ (30) 

The area of ICECUBE being A ^ 10^ m^, this gives Njce = l(tobs/year). That is, the number 
of neutrino detection events per year is approximately $^(1 TeV)/10^^ (see, e.g., Crocker, 
Meha, and Volkas 2000, 2002, 2005; Crocker, Fatuzzo, Jokipii et al. 2005). 

More precisely, the number of neutrino detection events in a particular device is 

//•year 
dE, / dt Area[E„ e{t)] $(E,) x Pdctect(^.) Attn[E„ e{t)] , (31) 

which includes attenuation of the flux as it passes through Earth, the effective muon area — 
the preceding depending on the location of the source defined in terms of its angle from 
the nadir 9 — and the probability of neutrino interaction and subsequent detection. For 
the ICECUBE detector, these quantities are specified in detail in the Preliminary Design 



Document available from the ICECUBE website (http://icecube.wisc.edu/ ). See also Ahrens 
et al. (2004). 

Carrying this out gives an event rate of 0.3 detection events per year from Coma via the 
ICECUBE detector. More optimistic are the detection rates should the 7-ray bright Abell 
clusters be confirmed: these range between 30 upcoming events per year from Abell 1914 in 
the most optimistic case of a Sp = 2.1 spectral index, to 8 x 10"^ events per year from an 
extremely strong shock in Abell 85 (for the most pessimistic case). 

Meanwhile, the prospect for using Auger to directly probe neutrino production via 
P7 collision processes in clusters is possible, but it would require very favorable system 
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parameters: even assuming the weakest possible shock (and the most optimistic duty cycle 
of 15%), yields only a few detection events per year, with Abell 1914 representing the most 
promising source. 

Table 2, summarizing these results, as calculated on the basis of the results presented 
by Miele et al. (2006), is given below; here, the fields separated by a semicolon indicate 
decreasing proton spectra, from Sp — 2.1 to 2.5, in steps of 0.1. At the ~EcV energy scale 
of the p7 neutrinos, other experiments may offer apertures superior to that of Auger. In 
principle, the ANITA experiment (Barwick et al. 2006), in particular, has a considerably 
larger neutrino aperture than Auger but, unfortunately, restrictions on the positions of po- 
tential sources (which must be located within a narrow band of declination [—9°, 19°] to be 
detectable) exclude all the potential cluster neutrino source considered here. 

4. Detections from Known Radio Clusters 

Let us suppose that 7-ray detections from Abell clusters 85, 1914, and 1758 are due to 
71^ decay. We may then use multiwavelength observations to set a very solid prediction for 
their neutrino flux, depending only on the spectral index of protons — and their maximum 
energy. These are shown in Fig. 4 for a range of proton spectral indices consistent with 
Fermi acceleration and the synchrotron observations. 

The spectral slope of radio emission for Abell clusters 85, 1914, and 1758 are > 1.5, 
1.19, and 1.13, respectively (Ensshn et al. 1998). Abell 85's value of a = 1.5 would require 
a proton spectral index above that which could be produced by shock acceleration. This 
could indicate either variability in measurements made years apart, or the influence of several 
bright radio galaxies within the cluster. Here we take the maximum value reasonable for 
shock acceleration, Sp = 2.5. The other values indicate proton injection spectra near Sp — 2.3, 
matching the cooled spectrum of Coma (see Eqs. 13-15). 

Of the three, Abell 85 appears to be the only cluster (Durret et al. 2005) for which 
temperature maps exist to create a radial gas profile, although Chandra detections of 1914 
do exist (Govoni 2005). For Abell 85 the core gas density is 7.7 x 10~^ cm~^, while in the 
other two clusters we assume it is the same as in Coma (ue = 3 x 10~^ cm~^). The flux 
of particles for each cluster scales as its relative redshift to Coma, while the active volume 
scales as the largest linear scale (LLS). In every case we take a nominal B = l/xG, for lack of 
better knowledge. This results in a strikingly similar value for the CR proton number density 
required to fit the radio flux (Fig. 3): in every case (including Coma) it is within a factor 
of two of ncR — 10"^ cm~^, where n(p) — ncRP"^"- This is very near energy equipartition 
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with the background thermal electrons. This balance in energy between the two species may 
be just a coincidence, but more likely it represents a dynamic coupling between the protons 
and the ambient electrons, such that an effective energy transfer occurs between the two (see 
Wolfe and Meha 2008). 

In each case, 7-rays from 7r°-decay dominate those from nonthermal bremsstrahlung by 
the synchrotron-emitting electrons. It is therefore impossible to describe both the Coma 
cluster and these 7-ray bright clusters within a simple secondary model, since Coma should 
then be a copious source of 7-rays. The situation is not improved if one assumes nonthermal 
bremsstrahlung from the electrons produces the 7-rays, since this would (via inverse Compton 
scattering) produce an X-ray excess not observed in Abell 1785 and 1914 (Abell 85, however, 
displays such an excess). 

De Marco et al. (2006) have advanced the notion that photopion production between 
cosmic-ray protons and infrared light emitted by galaxies within the cluster may dominate 
over proton-CMB interactions. We agree with this general result, but find some difference in 
the details. Again, pion production from p7 interactions is characterized by a hard cut-off 
at the threshold energy, followed by a spectral index one greater than the proton index, 
— Sp + 1. The three-body decay then gives an extra degree of freedom, allowing neutrinos 
to be produced at a constant emissivity until threshold. The De Marco et al. (2006) (their 
Fig. 5) and our own calculation include all of these features, except the location of threshold. 
The energy of the galactic IR emission under consideration (Lagache et al. 2003) is peaked 
a decade above the CMB peak; with the threshold for 7 -|- p — > p -|- tt being 

_ {rrip + m^f - ml 

^thres — , yO'^) 

the Wein peak at 1.2 x 10"^ eV gives a CMB threshold of 5.8 x lO"*^^ eV, while the IR peak 
at 1.03 X 10~^ eV gives a threshold of 6.6 x 10^^ eV (later line emissions are suppressed by 
the steepness of the proton distribution). Yet the threshold of the De Marco et al. (2006) 
calculation is ~ 10^^ eV. This shift in the emission peak means neutrino production begins 
at a point in the proton distribution some 100 eV in energy above that given in De Marco 
et al. (2006), which with n{p) ~ Ep^'^ corresponds to a factor of 10~^ for the neutrino fiux. 

Finally, we comment on the possibility of calculating diffuse cluster 7-ray and neutrino 
fluxes, using the Coma cluster as a template. In this case it is helpful to look at those 
radio halo/relic clusters which are known to exist, assuming that bright sources of 7-rays or 
neutrinos will have first been observed in one of the radio surveys. If cosmic-ray production 
is dominated by cluster mergers, the number of CR protons, and thus the radio luminosity 
(as well as neutrino and 7-ray production), should roughly scale as a cluster's mass, i.e., 
Lp ~ GM'^/Rt^eige ~ {M/Mcoma.y^^ (De Marco et al. 2006). As a given cluster's mass is 
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virially connected to its temperature, we may predict the flux of cosmic rays and neutrinos 
as a function of that cluster's relative temperature and distance with respect to Coma. In 
this case, M ~ 5T / G{k / ump) , where we take the virial radius as scaling with the largest 
linear size of a cluster. The M — T relation has also been derived from specific Chandra 
observations (Kotov & Vikhhnin 2006) as M = A{T/5f x lO^^M© with A = 1.7 and (3 = 1.5. 

Should the population of radio halo/relics be powered by gravitational infall at the 
cluster's current mass, we expect radio power to scale with the 5/3rds power of cluster 
mass — but it does not (Fig. 5). Using the Chandra M — T relation gives a correlation 
between (M/Mcoma)^^'^ and Ljadio of 0.46; the virial value improves this to 0.61. In no case is 
the relatively larger proton density of a particular cluster greater than (zcoma,/ z)"^ , meaning 
that none of the known radio clusters could contribute a neutrino or 7-ray flux above that 
of Coma's. The most concrete estimation for the diffuse flux — the weighted sum (by radio 
luminosity) of the 30 or so known radio clusters — is approximately 9/4ths that of Coma 
alone. Of these, 7-ray bright double clusters Abell 1758 and Abell 1914 (Fig. 5, square and 
diamond, respectively) are among the most dominant. The entire discussion is moot if the 
identification of Abell 85, 1758, and 1914 as 7-ray bright clusters is confirmed, in which case 
these three clusters alone dominate the diffuse flux despite their relatively low mass; in this 
case the importance of ongoing, violent mergers is emphasized. A generalization which would 
allow a sohd diffuse flux to be calculated appears to depend on a present lack of knowledge 
about the basic mechanisms of cluster halos. 



5. Conclusions 

Observations of both 7-rays (7pp) and neutrinos (both Upp and Up^,) in clusters provide 
answers to a series of yes-or-no questions. If both jpp and Upp are observed from a cluster, 
then proton-proton induced cascades are certainly at work and their contribution to the 
radio halo can be evaluated — as well as potential models for the creation of viable proton 
populations with total energies > 10^^ erg. If 7^^ is detectable by either EGRET or GLAST 
and Upp is not observed, then jpp is not likely to be caused by secondary decay (this is not 
a flrm 'no' as the source in question may occupy an unfavorable position in the night sky, 
or the proton population may experience a coincidental cutoff just after 100 MeV) and new 
models would need to be considered. An Auger Up^ detection would conflrm our expectation 
that very energetic cosmic-ray photons interact with the CMB to produce debris particles, 
thereby losing their energy. This would add some support to the view that cosmic ray events 
above 6 x 10^^ eV cannot be hadronic if the incident particles originate from beyond the 
GZK limit. Finally, if Upp is observed without Vp^, we may assume that the spectral index 
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of protons is steeper than 2.1, or that the cluster is incapable of accelerating protons to this 
energy. 

The combination of X-ray, radio, 7-ray, and neutrino observations on galaxy clusters 
would give a vital new perspective on the mechanisms producing high-energy particles in 
clusters. A simultaneous detection of 7-rays and neutrinos from pp scatterings in a clus- 
ter would confirm conclusively that secondary production is at work, and provide a new 
probe for cosmic-ray production in the MeV and TeV ranges, as well as shock acceleration in 
clusters and astrophysically in general. A detection (or upper limit) of p7 neutrinos, mean- 
while, would demonstrate whether clusters are capable of accelerating protons to ultra-high 
energies, and, indeed, how particle physics tested in colliders evolves towards these energies. 

A confirmation via GLAST of Abell 85, 1758, or 1914, as 7-ray bright would confirm 
them as strong neutrino sources, potentially detectable in the ICECUBE and, perhaps. Auger 
experiments. If a nonthermal proton population approaching energy equipartition is demon- 
strated, these then form a dominant source of pressure in clusters, affecting structure and 
evolution. Of course, the fact that neutrino detections are only possible near equipartition 
means that there may be only a narrow window of opportunity for actually measuring a 
neutrino flux, or maybe none at all. But we find it promising that an interesting predic- 
tion may be made even without requiring a super-equipartition proton population. And, 
anyway, a near-equipartition situation between the protons and the background electrons is 
suggested by a dynamic coupling between these two populations, as discussed in Wolfe and 
Melia (2008). Taken together with evidence of ongoing mergers, these observations would 
give new insight into the period of cluster formation, beginning some > 9 billion years ago, 
and what role mergers play in cosmic-ray production. 
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Table 1: Properties of Clusters Containing Halo and/or Relic Candidates 
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Table 2: Neutrino Detection Event Rates 



Name 
Abell 85 

Abell 1914 0.2; 0.8; 2; 8; 30 

Abell 1758 0.1; 0.4; 1; 5; 18 

Coma 0.3 



Auger (l/decade) 
0.08; 0.3; 1; 3; 13 7e-5; 9e-4; le-2; 2e-l; 2 

9c-5: le-3; 2e-3; 2e-l; 3 
5e-5; 7e-4; 9e-3; le-1; 2 
0.03 



IC ui)coniiiig (l/}'ear) IC downgoing (l/}'ear) 
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Fig. 1. — Cross-section for charged pion production via p7 scattering, and its fit. dashed- 
the fit for the A resonance; dotted-fit to the e' ~ 0.7 GeV resonance; solid-overall fit; stars 
data from Mucke et ah (2000). 
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Fig. 2. — Gamma rays (left panel) and associated muon neutrinos (right panel) from the 
Coma cluster. In the left panel, a flux of L^^^s = 0.3 is given by pion decay (thin solid), 
which wholly dominates nonthermal bremsstrahlung from synchrotron-emitting electrons 
(solid shaded). Also pictured are the detection limits for GLAST (1 year, dashed), EGRET 
(1 year, dash-dotted), VERITAS (thick solid), and MAGIC (5 a, 50 hrs, dotted). In the 
right panel we show Coma's associated neutrino flux for pp scattering (dash-dotted), p7 
scattering with the IRB (sohd shaded) and that with the CMB (thick solid). The 7-ray 
detection limits (dotted) and flux (thin solid) are shown for comparison. 
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Fig. 3. — Radio emission from tiie four clusters sets the scale for the number and index 
of cosmic rays. solid-Coma,; dashed-Ahell 1914; dotted-Ahell 1758; dash- dot- Ahell 85. In 
each of the four clusters examined here, the number of cosmic-ray protons is a few times 
10~^ cm~^, placing them just below energy equipartition with the thermal background. The 
implied similarity of the three clusters' environments (other than Coma) means that, in each 
case, 7-ray emission is dominated by 7r° decay. 
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Fig. 4. — A positive ~ GeV 7-ray detection directly implies a neutrino flux. In this case we 
display a model for Abell 1914 which describes its 7-ray flux as a secondary model (sohd, 
thin), leaving the precise spectral index {sp = 2.1 through 2.5, in steps of 0.1, moving down- 
wards) as an unknown. The EGRET and GLAST limits (dotted) are again reproduced as a 
guide. Neutrinos produced in concert with 7-ray emission from pp scattering (exclusively z/^ 
produced via muon decay, dash-dotted) dominate to several decades of TcV; those produced 
in p7 interactions with an infrared background (solid, shaded) peak at a few times 10^^ eV 
and dominate at a few times PeV; and p7 interactions with the CMB (thick solid) dominate 
above EeV. 
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Fig. 5. — Radio luminosity of radio halo/relics does not scale as their mass. It is therefore 
difficult to produce a template for calculating the total diffuse flux. The 7-ray bright radio 
halos Abell 1758 (box) and 1914 (diamond) are relatively more massive than other clusters, 
which may be why they are observed via EGRET. 



